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Final Project Report: 
Characterization and Design of 
Electromagnetic, Chemical and Thermal 
Transport Processes for Multi-Phase 
Systems, AFOSR Grant F49620-99-1-0009 

Robert P. Lipton. 


Abstract 

The work done in this project provides new mathematical tools 
necessary for the description of the qualitative behavior of composite 
materials. These tools are used to develop numerical algorithms for the 
design of composite structures. The first part of this project provides 
new analytical results necessary for understanding the macroscopic ef¬ 
fects of field concentrations arising from the composite microgeometry. 
These results form the basis of new a numerical method for the de¬ 
sign of elastic composites subject to stress constraints. The numerical 
methods provide rigorous guidance for the design of reinforced struc¬ 
tures for maximum specific stiffness and strength. Similar results are 
obtained for multi phase dielectric materials. A numerical method is 
developed to solve design problems requiring the control of the electric 
field in composite dielectrics. 

The second part of the project predicts novel physical phenomena 
for reinforced materials in the presence of imperfect bonding. These re¬ 
sults are unlike phenomena seen in more traditional treatments where 
the interface between constituents is treated as a ’’perfect bond”. The 
difference is that reinforcement size effects are predicted. New re¬ 
sults are found for imperfectly bonded nonlinear dielectric materials 
and reinforcement problems in linear elasticity. These results provide 
guidance for the robust design of particle and fiber reinforced elastic 
composite structures in the presence of imperfect bonding. 
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1 Research overview. 

During the last few years several technologically important applications 
have benefited from the use of functionally graded composite materials, 
see Glasser (1997) and Koizumi (1997). In many applications there is a 
separation of scales and the discrete entities forming up the microstructure 
exist on scales significantly smaller than the characteristic length scale of 
the loading. Under this hypothesis functionally graded materials are mod¬ 
eled using effective thermophysical properties that depend upon features of 
the underlying microgeometry. The effective thermophysical properties are 
given by effective constitutive laws relating average flux to average gradient, 
see Markworth, et. al. (1995). 

The primary problem of design of a functionally graded material is the 
determination of the optimal spatial dependence for the composition. This 
type of problem has generated much interest in the engineering community 
and is the topic of a rapidly developing literature, see Markworth, et. al. 
(1995), and Ootao, et. al. (2000). For many objective functions this type of 
problem has also received significant attention from both the applied math¬ 
ematics and structural optimization communities in the 1980s and 1990s 
under the headings of the homogenization method for topology optimiza¬ 
tion and structural optimization, see Bendose and Kikuchi (1988), Allaire 
and Kohn (1993), Lurie and Cherkaev (1986), Murat and Tatar (1985) and 
Cheng and Olhoff (1981). Homogenization methods applied to the design 
of composites for optimal structural performance can be found in the works 
of Sigmund and Torquato (1997), Fujii et. al. (2001) and Diaz and Lipton 
(2000). In all of these works the problem of determining the optimal spa¬ 
tial dependence for the composition is obtained through the use of effective 
constitutive relations. 

Motivated by the applications, this project treats the problem of op¬ 
timizing structural or dielectric properties subject to constraints on field 
quantities. These quantities include the local stress for elastic composites 
or the electric field for dielectric composites. In many applications it is of 
central importance to control fields inside composite structural components. 
Regions containing large field gradients are most likely the first to exhibit 
failure during service. For this type of problem the macroscopic effect of 
field concentrations that arise from the microstructure need to be accounted 
for. The concept of an effective constitutive law is by its self not sufficient 
to capture the effect of these concentrations. This requires new modeling 
beyond the notion of effective thermophysical properties. The objective of 
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this project is to point out a new class of macroscopic tensors relevant to the 
modeling of microscopic field concentrations and to present a methodology 
for the numerical design of functionally graded materials in the presence of 
constraints on field gradients. 

The second part of this project is concerned with modeling imperfectly 
bonded particle and fiber reinforced composites. During the previous pe¬ 
riod of AFOSR support this investigator developed a methodology to assess 
the effectiveness of reinforcement fibers or particles on the overall trans¬ 
port properties of composite materials. The approach has been successful 
in predicting the effect of particle size and shape on the enhancement of 
structural and thermal transport properties in the context of electric con¬ 
tact resistance, coupled heat and mass transport on the interface, highly 
conducting interfaces, and problems of torsional rigidity with imperfectly 
bonded fiber reinforcements. 

The work done here extends these techniques to imperfectly bonded 
nonlinear dielectrics, and to reinforcement problems in the context of two 
dimensional elastic systems. 

2 Grant sponsored projects completed. 

Homogenization of field fluctuations and stress constrained design 
of functionally graded engineering materials. 

The use of functionally graded materials is rapidly expanding and now 
includes structural, bio mechanical, and energy conversion applications. 
Functionally graded materials allow the designer to tailor the material’s 
microstructure in order to enhance structural performance. 

During the course of this project new homogenization results have been 
obtained that describe the effect of microscopic field concentrations on the 
macroscopic field intensity. It has been shown that these new quantities 
facilitate new numerical methods for the design of functionally graded ma¬ 
terials for maximum strength and stiffness. 

Functionally graded materials (FGMs) are characterized by microstruc¬ 
tures that are spatially variable on the macroscale. To fix ideas we consider 
a structural element made from a FGM consisting of elastic reinforcement 
fibers with stiffness tensor Ai embedded in a more compliant elastic matrix 
with stiffness tensor A 2 . The characteristic length scale of the microgeome¬ 
try relative to the size of the structure is denoted by e. In many applications 
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e is taken to be small, i.e., £ « 1. For a body load / and boundary traction 
g the stress o satisfies -diver = / and o n = g. The stress o is written as 
cr = o M + o™. Here o M is the macrostress and the macroscopic consti¬ 
tutive law is o M = C E e M , where e M is the average strain and C E is the 
effective elastic tensor. The macrostress satisfies —div a M = / and o^f = g. 
The microstress cr™ depends on e and captures the interaction between the 
microstructure and o M . 

We place a constraint on the stress o M + cr™ and for a prescribed load 
case, the design problem is to find the fiber distribution that renders the 
structural element the stiffest. We tackle the design problem for situations 
in which an extremely large number of fibers are used (e.g., thousands of 
fibers) and e « 1. Design problems of this sort are too large for numeri¬ 
cal solution using existing computational methods. Instead we look for an 
approximation to the problem that is computationally tractable and gives 
nearly optimal designs for the original problem. It has been shown by this 
investigator that such an approximation is found by taking the e = 0 limit 
of the original. However the correct identification of the e = 0 limit design 
problem is delicate. In fact when considering stress constrained problems 
the traditional micromechanical approach to the constitutive modeling of 
discreet systems fails to give the complete picture. In this limit, an accurate 
modeling of FGMs requires that the stress constraint include the effects of 
the microstress. The stress constraint can be local or of integral type and is 
written symbolically as C{o M + o™) < K. Recent work obtained during the 
course of this project shows that the stress constraint does not simply reduce 
to C(o M ) < K in the e = 0 limit. For a constraint of mean square type, 
i.e., / \o M + o™\ 2 dx < K, this investigator has shown that the homogenized 
constraint is given by 

lim J \o M + cr™\ 2 dx = J \o M \ 2 dx + J Qo M • o M dx < K (1) 

where Q is a new type of effective property and is given by 

Q(x) = |^(S E (x)A?V i C E (x)S E (A,x)| - I (2) 

Here S E (x) = (C E (x)) _1 is the effective compliance and V l C E gives the 
change of the effective elasticity due to a change in the stiffness of the i th 
component material. The convergence of the quantity |cr™| 2 is given by 

|a e ™| 2 -> Q(x)cr M (x) : cr M (x) (3) 
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in the sense of distributions and one has the estimate 

Q(x)cr M (x) : cr M (x) < liminf sup |cr™(x)| 2 (4) 

xinQ 

for each point x in the domain. It is evident from (1) that new effective 
properties beyond C E are required when considering the e = 0 limit of 
stress constrained problems. 

The design variables for the e = 0 design problem are obtained by ex¬ 
changing the original design variables with averaged design variables such 
as the local density of fibers. To illustrate the basic ideas we suppose that 
the fiber microstructure consists of a locally periodic arrangement of long 
parallel fibers with circular cross section. The fiber diameters are allowed 
to change over the macroscopic dimensions. For this case it has been shown 
that the design variable for the (homogenized) e = 0 design problem is the 
local density of fibers 9. The formulas for the effective tensors Q — Q(9) and 
C E = C E (0) are determined numerically through the solution of suitable 
microscopic problems. The resulting homogenized optimal design problem 
is one of maximization of the structural stiffness over 9 subject to the ho¬ 
mogenized stress constraint given by (1). The design problem is solved 
numerically using the method of steepest decent together with the finite ele¬ 
ment method. Numerical results have been obtained for the case of torsional 
loading for functionally graded fiber reinforced rods. 

The results reported in this project will appear in the following journals: 

1. R. Lipton, Design of functionally graded composite structures in the 
presence of stress constraints, International Journal of Solids and Struc¬ 
tures. To appear in 2002. 

2. R. Lipton, Relaxation through homogenization for optimal design prob¬ 
lems with gradient constraints, Journal of Optimization Theory and 
Applications. To appear in 2002. 

Characterization of electric-field fluctuations in random compos¬ 
ites. 

The theory of composite materials has for the most part focused on 
the characterization of effective transport properties that relate average flux 
fields to average gradient fields. In the context of heterogeneous dielectric 
materials this type of effective property is known as the effective dielectric 
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constant. The effective dielectric constant e e gives the linear relation be¬ 
tween the average electric field E and the average electric displacement D, 

he., ___ 

D = e e E. 

However, it should be pointed out, that a heterogeneous dielectric material 
exhibits a hierarchy of’’effective properties,” beyond s e . The first effective 
property in this hierarchy is the covariance tensor o. This effective property 
delivers the mean square fluctuation of the electric field about its average 
value E. To make this precise, the local electric field in the composite sample 
Cl is denoted by E(x). The local electric field is related to the local electric 
displacement D(x) by D(x) = e(x)E(x), where e(x) is the local dielectric 
constant and divD(x) = 0. Then cr is defined by 

< ’ 6 ' 6= M/ n |E(3l) - i:|2lix ' 

Higher order “effective properties,” are defined through higher moments 
of the local electric field fluctuation E(x) — E, i.e., 

i/ n |E(x)-ENx, P >2. 

It is evident that a and e e are correlated as they are generated from 
one and the same composite sample. For two phase composites made up 
of dielectric constants e\ and e 2 , it is shown that the relation between the 
effective dielectric constant and the covariance tensor is given by 

O = (e e - Sil)/ei - — — —d S2 e e (ei,£2)- 

Here we have written s e = £ e (£i,£2) to emphasize its dependence on the 
dielectric constants of the constituent materials. The second term on the 
right hand side is the derivative of the effective dielectric constant with 
respect to £ 2 - 

Higher order effective properties provide information on the variation of 
the local electric field that is not captured by the effective dielectric constant. 
The higher order effective properties reveal the presence of regions of high 
field intensity inside the composite. These regions are most often the first 
to suffer dielectric breakdown during service. 

Unfortunately, the higher order effective properties can not be obtained 
directly through simple boundary measurements. While the effective di¬ 
electric constant can be easily measured by subjecting a composite sample 
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to a uniform electric field and measuring the current passing through the 
boundary of the sample. In this project we have developed bounds on the 
covariance tensor of the electric field that are given in terms of the effective 
dielectric constant. This allows one to obtain estimates for field concentra¬ 
tions within the composite sample from simple boundary measurements. We 
exhibit microstructures for which these bounds are optimal. These bounds 
are used to recover bounds on the covariance tensor when only the phase 
volume fractions and the two-point correlation function are available. For 
isotropic composites we obtain a lower bound that is the most restrictive 
one in terms the volume fraction occupied by each phase. Lastly, we provide 
tight upper and lower bounds on the mean square field fluctuation when only 
the volume fractions of the two dielectric materials are known. Denoting the 
volume fraction of the £2 dielectric by 6 2 the bounds are given by, 


0 < crijEiEj < U(6 2 ,E), 


(5) 


where U{9 2 , E) depends upon the contrast A = £ 2/^1 and is given by 


U{0 2,E) 


(02f(h)) |E| 2 , if e 2 < 1 -h, 
(02/(1 - 02)) |E| 2 if 0 2 > 1 -h. 


Here h = 1/(A — 1) and the function f(z) is defined by 


m 


z 

(ih + z ) 2 ' 


These bounds are the best one can find on the covariance knowing only the 
volume fractions of the component phases. They are the analog of the well 
known harmonic mean arithmetic mean bounds for the effective conductiv¬ 
ity. Extremal sequences of configurations that attain the bounds are shown 
to be given by the well known finite rank laminate microstructures. The 
lower bound is attained by laminates of the first rank with layers oriented 
parallel to the applied field. On the other hand the upper bound is satu¬ 
rated by laminates of first or second rank depending on the magnitude of 
the contrast. The method for obtaining bounds developed in this project 
relies on the explicit computation of the convex hull of a curve. The results 
reported here have appeared in the following journals: 


1. R. Lipton, Optimal inequalities for gradients of solutions of elliptic 
equations occurring in two-phase heat conductors, SIAM J. Mathe¬ 
matical Analysis , 32, 1081-1093 (2001). 
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2. R. Lipton, Optimal bounds on field fluctuations for random compos¬ 
ites, J. Applied Physics , 88, 4287-4923 (2000). 

3. R. Lipton, Optimal bounds on field fluctuations for random compos¬ 
ites: Three dimensional Problems, J. Applied Physics , 89, 1371-1376 
( 2001 ). 

Optimal containment and control of dc electric fields and electric 
displacement fields. 

Over the last 25 years several methods have been developed for the 
design of composites with minimum energy dissipation for given source dis¬ 
tributions and boundary conditions. These design problems appear in many 
contexts and include the use of dielectric and elastic composites. However, 
many practical problems require that field quantities be controlled directly. 
Unfortunately none of the computational strategies developed for minimum 
energy dissipation can be used. The basic reason is that the objective func¬ 
tions appearing in problems where the fields need to be controlled directly 
are not continuous with respect to G-convergence. Recently A. Velo and 
myself have been able to circumvent this obstacle by computing the G- 
continuous extension of a class functionals used in the direct control of field 
quantities through the use of layered microstructures. 

Based upon this extension a computational methodology has been devel¬ 
oped for the design of functionally graded composites used for the contain¬ 
ment and control of dc fields. We start with a two-phase dielectric material 
occupying the design domain fh A resource constraint is placed on the 
amount of the better dielectric that can be used in any admissible design. It 
is assumed that a prescribed charge density has been distributed inside the 
design domain. It is desired that in some sub-region D of Q, that the electric 
field have an intensity and direction specified by a prescribed target field E. 
The goal is to find the configuration of the two dielectrics so that the actual 
electric field is as close as possible, in the mean square norm, to the target 
electric field. In this project we have introduced a tractable method for the 
numerical computation of minimizing sequences of configurations of the two 
dielectrics. The configurations are associated with materials with graded di¬ 
electric properties that may exhibit a fine scale structure composed of layers 
of the two dielectrics. The computational results provide rules of thumb for 
the design of graded dielectric materials. These results are reported in: 

1. R. Lipton and A. Velo, Optimal design of gradient fields with appli- 
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cations to electrostatics, in Nonlinear Partial Differential Equations 
and Their Applications, College de France Seminare Series in Applied 
Mathematics, edited by P.G. Ciarlet et P.-L. Lions. Elsevier-Gauthier- 
Villars , to appear in 2002 . 

The effect of the interface on the dc transport properties of non¬ 
linear composite materials. 

Nonlinear dc electric conductivity is a property intrinsic to many ce¬ 
ramic materials. More generally, nonlinear inhomogeneous materials are 
pervasive, appearing in applications ranging from transient voltage protec¬ 
tion to the selective absorption of solar energy. This work investigates the 
effect of the interface on the overall dc electric properties of nonlinear com¬ 
posite conductors. A composite consisting of a monodisperse suspension of 
spheres embedded in a matrix is considered. It is supposed that the interface 
separating the spheres and matrix is highly conducting. A critical applied 
voltage is found for which the electric potential inside the sample is the same 
as for a sample of identical shape containing no spheres whatsoever. At the 
critical voltage, the overall electric current passing through the sample is the 
same as in a homogeneous conductor. On the other hand, in the presence of 
an electric contact resistance at the interface, we show that there is a critical 
applied dc current density for which the current density inside the sample 
is the same as for a sample containing no spheres. The overall electric field 
in the sample corresponds to that associated with a homogeneous conduc¬ 
tor made from matrix material. These effects are shown to be independent 
of the location of the spheres within the sample. Moreover, this effect is 
independent of the concentration of spheres in the sample even beyond the 
onset of interface percolation. 

It is demonstrated from first principles that this phenomenon occurs for 
a wide range of nonlinear constitutive behavior. Indeed, it is shown that 
this phenomenon can occur when the potential energy density of each phase 
is a convex function of the magnitude of the electric field. This requirement 
naturally includes local constitutive relations associated with nonlinear be¬ 
havior of the form 

j = 7 |E|* E 

where j is the local current density and E is the local electric field. Here the 
nonlinear susceptibility 7 takes different values inside each phase. 

This work is joint work with David Talbot and has appeared in 
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1 . R. Lipton and D.R.S. Talbot Journal of Applied Physics, 86, August 
1 , 1999, pp. 1480-1487. 

Effect of interfacial bonding on fiber reinforced shafts subject to 
antiplane shear. 

Fiber reinforced materials are often the materials of choice for structural 
components appearing in aerospace and infrastructure applications. The 
fibers can be perfectly bonded to the matrix or bonded to the matrix with a 
bond that is stiffer than either the fiber or matrix phase. Over time several of 
the fibers may become imperfectly bonded to the matrix. In this project we 
consider a shaft reinforced with long fibers subjected to antiplane shear. We 
outline conditions for which the effects of the interface overcome the elastic 
properties of the fiber reinforcement. We apply a new geometric criteria that 
indicates when the imperfect bond compromises the elastic properties of the 
reinforcement fiber. The criteria is given in terms of the bond stiffness and 
the surface traction to bulk stress quotient of the fiber cross section. This 
criteria was introduced by this investigator during the previous AFOSR 
funding period. This criteria is used to provide new rigorous rules of thumb 
for selecting the fiber size so that the stiffness of the structural component 
is greater than the unreinforced matrix. 

On the other hand in the applications it is anticipated that one only has 
statistical data on the type of bonds and bond stiffness associated with a 
fiber reinforced structure. We present new mathematically rigorous guide¬ 
lines for the design of the overall stiffness of fiber reinforced shafts when 
only partial statistical information is available. These results are applied 
to obtain new size effects for distributions of fibers with circular cross sec¬ 
tions. These results are consequences of new upper and lower bounds on the 
effective compliance tensor associated with antiplane shear loading. 

These results are applied to the design of a steel reinforced concrete 
column subject to antiplane shear. Given that a certain percentage of the 
reinforcement rods are imperfectly bonded we are able to predict a range of 
rod diameters for which no increase in the overall shear stiffness is possible. 
The utility of the methods presented in this project lie in the fact that 
simple but very general design rules are obtained without the use of intensive 
computation. Naturally more detailed relationships between fiber geometry 
and overall structural properties will require significant numerical effort. 

This work has appeared in the following journals. 

1. R. Lipton, Effect of interfacial bonding on fiber reinforced shafts sub- 
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ject to anti plane shear, International Journal of Solids and Structures , 
38, 369-387 ( 2001 ). 

2. R. Lipton, Effect of interfacial bonding on fiber reinforced materials. 
In the Proceedings of the 13th ASCE Engineering Mechanics Division 
Specialty Conference, Johns Hopkins University, Baltimore, June 13- 
16 1999. (CD-ROM) 

Reinforcement of elastic structures in the presence of imperfect 
bonding. 

In this work we consider the two dimensional problem of plane strain 
in the presence of imperfectly bonded elastic reinforcements. The approach 
taken here extends the methodology developed by this investigator for scalar 
problems to treat reinforcement problems described by the system of linear 
elasticity. A new geometric criterion on the shape and size of the elastic 
reinforcement is found that determines when the effects of imperfect bonding 
overcome the benefits of the reinforcement. The criterion is given in terms of 
a new type of eigenvalue problem posed on the surface of the reinforcement. 

Both the matrix and reinforcement are assumed to be made from isotropic 
elastic materials. The Lame constants for the reinforcement and matrix 
phases are given by /i r , A r and X m respectively. Here the reinforcement 
is stiffer than the matrix. This condition is equivalent to the requirement 
that k t > K m and \i r > ^u m , where K r and K m are the plane strain bulk 
moduli of the reinforcement and matrix respectively, i.e., n r — (i r -V \ r and 
K>m — l^m X rn . 

On the other hand the interface separating the reinforcement from the 
matrix is imperfect. This is associated with a discontinuity in the displace¬ 
ment across the interface. The constitutive model for the interface is taken 
to be of the spring layer type. This model represents a first order approxima¬ 
tion to more general nonlinear interface models, e.g., the yielding interface 
models or shear lag models. For the spring layer model the jump in dis¬ 
placements is proportional to the traction at the interface. The constant 
of proportionality is denoted by /3 and has dimensions of stiffness per unit 
length. 

We define the relative compliance 7 of the interface and reinforcement 
system by 

r 1 _ r 1 

(2m™)- 1 - (2Mr)- 1 ’ (2 k™)- 1 - (2/Cr)- 1 




7 = mm 


( 6 ) 
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where 7 has dimensions of length. 

In this work we show that the interplay between the imperfect interface 
and the effect of the stiff reinforcement is mediated through a geometric 
quantity o related to the shape and size of the reinforcement. We suppose 
the structure is reinforced by several domains containing stiff material. We 
focus our attention on one of the reinforcement domains S and examine its 
effect on the overall stiffness of the structure. The quantity o associated with 
the reinforcement £ is written as cr(£). The quantity ct(£) is the stationary 
value of the Rayleigh quotient: 


a(£) 

where, 

C 


min 

<p£C 


IdE dstyrt dstyrt dl 


(7) 


|tp G tf 5/2 (£) 





Here d s denotes tangential differentiation on the curve <9E and A 2 is the 
biharmonic operator. The notation < p„ represents differentiation, i.e., <p,i = 
d Xi <p and tp^j = d x . x .(p, The parameter a has dimensions of inverse length. 
Conditions of stationarity for (7) deliver the eigenvalue problem 


A 2 0 = 

0, on £, 

(8) 

nidg<p,i = 

— aM n (<j >), on d£, 

(9) 

d s {tid 2 s (j),i) = 

-o{d s M s (4>) + Q{4>)}, on <9£. 

(10) 


Here is the bending moment of the reinforcement E given by M n (<f>) = 

riiTij(f) : ij and d s M s ((f>) + Q{<f>) is the Kirchoff shear force where M s (</>) = 
U^j<f>iij an d Q{4>) = d n Acj> = 0 . We have denoted normal differentiation by 
d n and the Laplace operator by A. Prom its definition we see that a is the 
largest constant C for which the inequality 

/ d s ip,i d s <p,i dl>C (p,ij (p,ij dx. ( 11 ) 


holds for all ip in the space C. For a disk of radius a one has that 



( 12 ) 


The overall stiffness of the reinforced structural cross section fi is in¬ 
versely proportional to the compliance energy of the structure associated 
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with a prescribed traction g. The criterion for when the imperfect interface 
overcomes the stiffening effect of the reinforcement domain S is given by 

Compliance energy inequality. 

If a satisfies, 

7 _1 < c( s ), (13) 

then for every load case p, the reinforcement does not reduce the compliance 
energy. We emphasize that this result is independent of the location and 
geometry of the other components of the reinforcement phase and applies to 
every load case g. 

The compliance inequality naturally implies a reinforcement size effect. 
Here the length scale is set by 7 . The effect of reinforcement size is seen 
clearly when E is a disk of radius a. We have the following size effect for a 
circular reinforcement. 

Size effect for a circular reinforcement. 

If the reinforcement is a disk of radius a, and 

a < | t ) ( 14 ) 

then for every load case g , the reinforcement does not reduce the compliance 
energy. 

This result gives a rigorous rule of thumb for the selection of the size 
of a reinforcement disk, namely: Only disks of radius greater than ( 2 / 3)7 
can provide reinforcement. This statement holds true independently of where 
the disks are placed in the structure. This result is in striking contrast to 
what is seen when there is perfect bonding between structural materials. For 
that situation, the addition of an infinitesimally small disk of stiffer material 
always reduces the compliance energy. 

These results are directly applicable to problems of optimal compliance 
design for plane strain problems. A prototypical problem is the optimal 
design of a structure reinforced with disks of different radii. Each disk has 
compliance C~ l and the matrix has compliance C” 1 . Here we suppose that 
each disk is made from stiffer material, i.e., C~ l < C^ 1 ■ The class of 
admissible designs is given by the set of all reinforcements consisting of a 
finite number of non-intersecting disks. We restrict the joint area of the 
disks to be less than a prescribed area fraction 0 r of the design domain ft. 
However, no lower bound is placed on the size of the disks nor do we place 
a constraint on the number of disks appearing in any design. We show in 
this project that all energy minimizing configurations of disks can be found 
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among those that contain disks of radii greater than or equal to 57 or no 
disks at all. It is evident that minimizing sequences made from progressively 
finer suspensions of disks (i.e., homogenized designs ) can be excluded. 

This work has appeared in 

1. R. Lipton, Reinforcement of elastic structures in the presence of im¬ 
perfect bonding, Quarterly Journal of Applied Mathematics (2001), 
pp. 353-364. 

3 Work at Air Force laboratories. 

I have been interacting with scientists in the Nonmetallic Materials Division 
at Wright Patterson Airforce Base since May of 2000. During the Summer of 
2001 I met with scientists at the Air Force Research Laboratory at Wright 
Patterson Air Force Base and it was agreed that the design of function¬ 
ally graded materials fits well with the research interests and goals of the 
Nonmetallic Materials and Ceramics Divisions. Much of the basic results 
on functionally graded materials obtained during the course of this project 
will provide the foundation for my interaction with the Air Force Labora¬ 
tory. Future projects related to the design of composite sandwich structures 
and reinforced ceramics will be done in collaboration with Dr. Iarve of the 
Nonmetallic Materials Division and Dr. Kerans and Dr. Parthasarathy of 
the Ceramics Division. A third project motivated by results obtained un¬ 
der this grant will be to compute bounds on elastic field fluctuations for 
realizable classes of random geometries that have been developed by Dr. 
Jonathan Spowart of the Metallic Composites Division at the Air Force 
Research Laboratory. 

4 Publications from grant sponsored activity given 
a Featured Review in Mathematical Reviews 

Optimal fiber configurations for maximum torsional rigidity, Archive for 
Rational Mechanics and Analysis, 144 (1998), pp. 79-106. 

Mathematical Reviews reference number: 2000i:74075. 
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5 Publications resulting from grant sponsored ac¬ 
tivity. 

1 . Optimal inequalities for gradients of solutions of elliptic equations oc¬ 
curring in two-phase heat conductors, SIAM J. Mathematical Analysis , 
32, 1081-1093 (2001). 

2 . Optimal bounds on field fluctuations for random composites, J. Ap¬ 
plied Physics , 88 ? 4287-4923 (2000). 

3. Optimal bounds on field fluctuations for random composites: Three 
dimensional Problems, J . Applied Physics , 89, 1371-1376 (2001). 

4. Effect of interfacial bonding on fiber reinforced shafts subject to anti 
plane shear, International Journal of Solids and Structures , 38> 369— 
387 (2001). 

5. Reinforcement of elastic structures in the presence of imperfect interfa¬ 
cial bonding, Journal of Engineering Mechanics , 127 . 667-671 (2001). 

6 . Design of functionally graded composite structures in the presence of 
stress constraints, International Journal of Solids and Structures . To 
appear in 2002. 

7. Relaxation through homogenization for optimal design problems with 
gradient constraints, Journal of Optimization Theory and Applica¬ 
tions. To appear in 2002. 

8 . Optimal design of gradient fields with applications to electrostatics, 
in Nonlinear Partial Differential Equations and Their Applications, 
College de France Seminare Series in Applied Mathematics, edited by 
P.G. Ciarlet et P.-L. Lions. Elsevier-Gauthier-Villars , to appear in 
2002. (Ani Velo - coauthor.) 

9. Optimal material layout in three-dimensional elastic structures sub¬ 
jected to multiple loads, Mechanics of Structures and Machines , 28 
(2000), no. (2&3), 221-238. (Alejandro Diaz - coauthor.) 

10. Bounds for the effective conductivity of a composite with an imper¬ 
fect interface, Proceedings of the Royal Society of London, 457 (2001), 
1501-1547. (David Talbot - coauthor.) 
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11 . An isoperimetric inequality for the torsional rigidity of imperfectly 
bonded fiber reinforced cylinders, Journal of Elasticity , 55, (1999), 

pp. 1 - 10. 

12. Variational methods, bounds and size effects for two-phase composites 
with coupled heat and mass transport processes at the two phase in¬ 
terface, Journal of the mechanics and Physics of Solids, 47 (1999), pp. 
1699-1736. 

13. The effect of the interface on the dc transport properties of nonlin¬ 
ear composites, Journal of Aplied Physics, M (1999), pp. 1480-1487. 
(David Talbot - coauthor). 

14. Reinforcement of elastic structures in the presence of imperfect bond¬ 
ing, Quarterly Journal of Applied Mathematics (2001), pp. 353—364. 

15. Effect of interfacial bonding on fiber reinforced materials. In the Pro¬ 
ceedings of the 13th ASCE Engineering Mechanics Division Specialty 
Conference, Johns Hopkins University, Baltimore, June 13-16 1999. 
(CD-ROM) 
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